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I. Experimental Observation 



A. Fixed Voltage - Switch Closed (v = V ) 




Switch 




+ 



+ 



V 



A 



£ 



v 



As an insulating material enters a free-space capacitor at constant voltage 
more charge flows onto the electrodes; i.e. as x increases, i increases. 

B. Fixed Charge - Switch open (i = 0) 

As an insulating material enters a free space capacitor at constant charge, 
the voltage decreases; i.e. as x increases, v decreases. 

II. Dipole Model of Polarization 

A. Polarization Vector P = Np = Nqd (p = q d dipole moment) 



N dipoles/Volume (P is dipole density) 




6.641, Electromagnetic Fields, Forces, and Motion 
Prof. Markus Zahn 



Lecture 7 
Page 1 of 27 



Nc field 



Electric field E 



! © ! 



Electronic polariiatian 



F = ,,E 



Torqufl - d x?E 



Orientitiwi and ionic polarization 
lb) 



Figu re S- 1 A n electric d ipo!e con sisis of two charges of equal magnitude but op posite 
sign, separated by a small vector distance d («) Electronic polarization anws when the 
average motion of the electron cloud about its nucleus i* slightly displaced, (ft) Orien- 
tation polarization arises when an asymmetric polar molecule tend* to line up with an 
applied electric field. If the spacing d aUo changes, the molecule has ionic polarization. 




d-da < 

Volume V 
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Qinsidev = -^qNd-da = }p P dV 

v y 

paired charge or 
equivalently 
polarization 
charge density 



Qinsidev = ~j> p • da = -J v • P dV = | Pp dV (Divergence Theorem) 

S V V 

P = qNd 
V-P = -p P 



B. Gauss' Law 

V-(s E) = P total = p u +p p = p u - V.P 



unpaired charge 
density; also 
called free charge 
density 



V.(s E + p) = Pu 

D = s E + P Displacement Flux Density 

V-D = p u 
C. Boundary Conditions 









I Ea, Pa, D a 






a 










f E t , P bf D b 


b 
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V • D = p u => j> D • da = | p u dV => n • [D a - D b ] = a su 

S V 

V«P = -p p => j>P-da = -Jp P dV => n.[p a - Pb] = -a sp 

s v 

v *(e E) = p u + p P => <j>e E-da = J (p u + p P )dV ^> n-8 [E a - E b ] = a su + a s 

S V 

D. Polarization Current Density 



AQ = qNdV = qNd-da = P-da 



[Amount of Charge passing through 
surface area element da] 



5AQ dP -7- 
di = — — = — -da 
p at at 



[Current passing through surface 
area element da] 



= J P «da 

\ 



polarization current density 



J P = 



ap 
at 



Ampere's law: 



V x H = J u + J P + 8 



a§ 
at 



^ ap aE 

= Ju + +8. 

at at 



= J U+ !(s E + F>) 



= Ju 



a_D 
at 



6.641, Electromagnetic Fields, Forces, and Motion 
Prof. Markus Zahn 



Lecture 7 
Page 4 of 27 



III. Equipotential Sphere in a Uniform Electric Field 




Mm 0> (r, 0) = -E D r cos 
<D(r = R,0) = 



[O = -E z = -E o rcos0] 



O(r,0) = -E o 



cos 



This solution is composed of the superposition of a uniform electric field 
plus the field due to a point electric dipole at the center of the sphere: 



dipole 



p COS 

47i8„r 2 



with p = 47te E R 3 



This dipole is due to the surface charge distribution on the sphere. 



a s (r = R, 0) = 8 E r (r = R, 0) = -s 



~dr 





\ 2R 3 






1 + — 3- 




r=R 


r 


r=R _ 


= 3s E 


cos 8 





cos 
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IV. Artificial Dielectric 



,_ v _ 8 v 

E= d' ^ =8E = ^ 



8 A 

q = ° s a = — v 



r _ q _ 8 A 
v~ d 



+ 









lLj; . 




•K.V7: 



(a) 



T 2r: 
d 

I 



(b) 



Figure fi.fi. 1 (a) Plane parallel 
capacitor with region between 
electrodes occupied by a dielectric, 
(b) Artificial dielectric composed of 
cubic array of perfectly conducting 
spheres having radius R and spacing 



Courtesy of Hermann A. Haus and James R. Melcher. Used with permission. 

For spherical array of non-interacting spheres (s >> R) 

P = 4 n s R 3 E Tz => P z = N p z = 4n s o R 3 E N 



N = l/s 

/s 



P = e„ 



4 71 



E = v|/ e 8 E 



v|/ e = 4 xc 



v 



V v ' 



y (electric susceptibility) 



D = s o E + P = e o [l + Ve ]E = eE 

s r (relative dielectric constant) 



8 = S r 8 = 8 [1 + Ve ]= 6 1 + 4JI 



3\ 
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V. Demonstration: Artificial Dielectric 



(£)- ampi 
osc. 



■if- 



\ I 1 




Figure 6.6.2 from the microscopic 
point of view, the increase in 
capacitance results because the 
di poles adjacent to the electrode 
induce image charges on the 
electrode in addition to those from 
the unpaired charges on the opposite 
electrode. 



.balanced 
amplifier 




step-up -p -f 

cage ~ 

Figure 6.6-3 Demonstration in which change in capacitance is used to measure the 
equivalent dielectric constant of an artificial dielectric. 



Courtesy of Hermann A. Haus and James R. Melcher. Used with permission. 
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Figure 6.6.4 Balanced amplifiers 
of oscilloscope, balancing 
capacitors, and demonstration 
capacitor shown in Figure 6.6.4 
comprise the elements in the bridge 
circuit. The driving voltage comes 
from the transformer, while v G is 
the oscilloscope voltage. 



Courtesy of Hermann A. Haus and James R. Melcher. Used with permission. 



. 8 A _ q s A 

|Ai| = co AC |V| = ^2. 
II II R, 

R=1.87 cm, s=8 cm, A= (0.4) 2 m 2 , d = 0.15m 
oo=2ti(250 Hz), R s =100 kn , V=566 volts peak 
AC=1.5 pf 

v = coACR s |V| 

= (2ti) (250) (1.5 x 10" 12 ) (10 5 ) 566 = 0.135 volts peak 
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VI. Plasma Conduction Model (Classical) 



dv + = - vp + 

m =qE-mvv — 

at n + 

dv_ E - Vp 
m =-qE-mvv 

" dt n_ 

p + = n + kT , p = n kT 

k=1.38 x 10" 23 joules/°K Boltzmann Constant 

A. London Model of Superconductivity [T -> 0,v ± -> 0] 

dv+ = dv- _ 

m + — =q + E , m_ — = -q_E 

J + = q + n + v + 3 = -q_ n v_ 

d3+ d / - x dv + (qJ) q + 2 n + 

= — q n v = q n = q n = — — - E 

dt dt^+ + + / 4+ + dt + + m + lJti^ 

dJ_ d / - x dv_ q 2 n — 

-dF = -dt( q - nv -) = - qn -^r = - qn - V^=^7 

CO 2 8 

co? = q+ n+ , co? = — ^ (co p = plasma frequency) 

p+ m + e p - m e 

For electrons: q.= 1.6 x 10" 19 Coulombs, m.=9.1 x 10" 31 kg 
n.=10 20 /m 3 , 8 = 8 * 8.854 xlfr 12 farads/m 



co p = 



q_ 2 n_ 



5.6x10" rad/s 



t = 



= ffi gxl0 10 Hz 

271 
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B. Drift-Diffusion Conduction [Neglect inertia] 



dy? ^ 

m+ w = q+ E ~ m + V + V + 



V(n + kT) - q - kT n 
— — '- => v = ^+ E Vn, 



m m,v,n 



dvi 

m — V— = -q E - m v v 
it 



V(n kT) - — q — k T 

_ i ^=>v = — — E Vn 

n m v m v n 



m + v + m +v+ 



- q 2 n - q kT „ 

J- = -q n v = — — - E + — Vn 

m v m v 

P + = d + n + , P_ = -q_ n_ 

J + = P+ n + E-D + V P+ 

J- = -p_ n_ E-D_Vp_ 



m 1 v i 



D = 



kT 

rn v. 



q_ n kT 

m v m v 



Y Y 

charge mobilities Molecular 

Diffusion 
Coefficients 



D + _ D _ kT 

n_ q 



= thermal voltage (25 mV @ T*300° K) 



Einstein's Relation 
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C. Drift-Diffusion Conduction Equilibrium (j+ = j_ = o) 

J + = = P+ n + E-D + v P+ = - P+ u + vo-D + v P+ 

J- = = -p_ n_ E - D_Vp_ = p_ \x_ VQ> - D_ Vp_ 



D -k T 

V<£ = — Vp + = V(lnp + ) 

p + n + q 



D k T 

V<£ = — — Vp = V(lnp ) 

P_ q 



P + = P e- qO/kT 



P. = -P„e + - /CT 



Boltzmann Distributions 



p (<d = o) = -p_ (<£ = 0) = p [Potential is zero when system is charge neutral] 



V 2 (D = ZP_ = ( p+ + p -) = zPjLr e -q*/kT _ e+q o./kTi - 2 P° C jnh qcl) 
8 8 8 L -I s kT 

(Poisson-Boltzmann Equation) 



Small Potential Approximation: -j^p- « 1 



. , qO q<B 
smh — — « - — 
kT kT 



V 2 $-^$ = 
skT 



V 2 <X> - = ; L d = Debye Length 

L d 2 d V2p q 
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D. Case Studies 

1. Planar Sheet 



<I>(x = 0) = V 



d 2 <t <t 



--^ = => <D = A x e x/Ld + A 2 e~ x/Ld 



dx 2 L 



B.C. <D(x -> +00) = 
O(x = 0) = V D 



O(x) 



" V e" x/L - x > 
"1 V e +x/Ld x < 
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r V q e -x/L 



do 

E ' = -*T = < 



d x > 



-V„ 



— s.e x/Ld x<0 



r _ 



dE x 



eV, 



L d 2 



e - x / L d 



EV 



L d 2 



g+x/L d 



x > 



x < 



2sV„ 



a s (x = 0) = e[E x (x = + )-E x (x = 0_)] = ^ 



2. Point Charge (Debye Shielding) 



0> 



V 2 <D r = 

^ y ^ 



1ST, 5<t^ 



r 2 dr 



1 d 2 



dr 



(r<X>) 



r sr 

E. Ohmic Conduction 



dr 2 1 7 L 2 







ro = a^"^ + /^e +r/Ld 
<D(r) = — ^— e" r/Ld 



J+ = P+ u + E-D + V P+ 
J = -p_ u_ E-D_Vp_ 

If charge density gradients small, then Vp + negligible => p + 

J = J + + J_ = ( p+ - p _ u_)E = p (u. + + u_)E = aE 



= "P. = Po 



o = ohmic conductivity 



J = a E (Ohm's Law) 



6.641, Electromagnetic Fields, Forces, and Motion 
Prof. Markus Zahn 



Lecture 7 
Page 13 of 27 



F. pn Junction Diode 



P= Na 1 " 


+ 
+ 


n = Nd 


n = n i ! " 

N fl ; " 


+ 
+ 


! p _ ni 2 



— I — h- 

-Xp 



-> x 



Xp 



qN A 



X n 



> X 

qN D x n - qN A Xp 
(Charge neutrality) 



Ex 
A 



-Xp 




-qN A (x+Xp) 

E 



E>! _ +qNp(x-x n ) 



C|lp = 



+qN A (x+Xp) 2 +<p p 
2e 



o = - 



qN D (x-X n ) Z +d'n 




On 



"Xp 



> x 



>j d> r - d) P =JSL | n 



L n' J 



In 



Hi 



n - N D =nie* 1 * ta ' CT 
^ q 111 ni 



x p = 



'2e\ Nd 



(o n -a>p) 



(N A +N D ) 



Xn _ 



-12 



(N A +N D ) 
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kT, IM. N n 

AO = (D - (D = In -j 

p q n, 2 



p 28 28 



= ^^(x n+ x p )= qN ° X " 2 fl3 
2e v n ~>c I m 



28 
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VII. Relationship Between Resistance and Capacitance In Uniform Media Described 
by e and a . 



Perfectly conducting 




Mathematical \ 
surface S 




Check: 



Parallel Plate Electrodes: R = — , c = — ^ RC = 8 / 

a A I /<* 
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Coaxial 
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VIII. Change Relaxation in Uniform Conductors 



V-J u+ ^=0 

at 



V • E = Pu / 



Ju = ctE 



a V • E + 
V 



3p u _ 



St 



= => 



^ + -P U =0 
St e Fu 



x e = s/cj = dielectric relaxation time 



St 



^ =0 



p u = p (F, t = 0) e-^ 



IX. Demonstration 7.7.1 - Relaxation of Charge on Particle in Ohmic Conductor 



A' 




Figure 7.7.1 Within a material having 
uniform conductivity and permittivity, 
Initially there is a uniform charge density 

in a spherical region, having radius 
a. In the surrounding region the charge 
density is given to be initially zero and 
found to be always zero. Within the 
spherical region, the charge density is found 
to decay exponentially while retaining its 
uniform distribution. 



Courtesy of Hermann A. Haus and James R. Melcher. Used with permission. 




Figure 7.7.2 The region between 
plane parallel electrodes is filled by 
a semi-insulating liquid. With the 
application of a constant potential 
difference, a metal particle resting 
on the lower plate makes upward 
excursions Into the fluid. [See 
footnote 1.] 



Courtesy of Hermann A. Haus and James R. Melcher. Used with permission. 
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Mathematical 






Surface S 

i 






Y 


- — ■ 




1 + > 


+ 


\ 




V 


I 


i A 






V 






\ + 




/ 




/ 











^ + ^=0^q = q(t = 0)e^ ( t . = s/) 
Partially Uniformly Charged Sphere 




Figure 3-21 An initial volume charge distribution within an Ohmic conductor decays 
exponentially towards zero with relaxation time r - e\a and appears as a surface 
charge at an interface of discontinuity. Initially uncharged regions are always un- 
charged with the charge transported through by the current. 
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Po 



Pu (t=0)=- 



r<R! 



r>R 1 



Q T = 3^ PO 



Pu (t) = r Po e 

■« 

_ 



r<R x 
r>R, 



(t. = e/o) 



Pore" 1 '" _ Qre-^ 



E r(r,t)= J 



3 8 4nsR, 3 

Q e-V'- 



< r < R x 



4 7c s r 



I _Q_ 



4 re e r 



Ri < r < R 2 
r > R, 



<t su (r = R 2 ) = e E r (r = R 2+ ) - e E r (r = R 2 _) 
_ Q 



= _y_^l_e-*.) 



X. Self-Excited Water Dynamos 

A. DC High Voltage Generation (Self-Excited) 




+ Q)Ring 



n drops/sec 



As the drops form near the rings. Image charges are Induced on the surface of the 
water, when the drops are completely formed (and Insulated firom one another) they 
carry a net charge, the sign of which Is determined by the sign of the charge on the 
respective rings. 
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Water drops Tall into ihc cans through cross-con netted wire loops. A potential 
difference of more than 20 kV between cans is spontaneously generated by the 
motion of the drops. For optimum operation the drops should form nearer to the 
rings than shown . This is accomplished by increasing the flow rate. 

Courtesy of Herbert Woodson and James Melcher. Used with permission. Woodson, Herbert H, and James R. Melcher. 

Electromechanical Dynamics, Part 2: Fields, Forces, and Motion. Malabar, FL: Kreiger Publishing Company, 1968. ISBN: 9780894644597. 
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nC iVl =C^ 



nC l \l 1 = C s V 2 



v 2 = V 2 e st 



n q V 2 = C s y 1 



nq 

Cs 



Y 

Det = 



nQ 

Cs 

j 



= 



nCY . nC 
— '- = 1 => s = ± — - 
Cs C 



8 root blows up 



Any perturbation grows exponentially with time 



B. AC High Voltage Self - Excited Generation 



Q>1 4= =r ©^ 



-nq v t 



-nq v 2 



-nq v 3 



c dv ^ 



dt 



c dv 3 



= C 



dt 
d^ 
dt 



v 2 = V 2 e s 



v 3 = V 3 e st 



nC Cs 



Cs 



nq Cs 

o nq 

J 



Vi 



V, 



= 



Y 

det = 
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(nQ) 3 + (Cs) 3 =0^s = f^|(-l)* 
Sj = -nQ/C (exponentially decaying solution) 



nq 

2C 



(-1) 



1/3 



1 

' 2 



[l ± V3 jj (blows up exponentially because s rea i >0 ; but also 
oscillates at frequency Sim ag * 0) 



XI. Conservation of Charge Boundary Condition 



i 




(a) 




TT 




(b) ■'' '=> 




St 




s ai - V 









n-pa -Jb] + ^CT SU = 
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XII. Maxwell's Capacitor 



Area A 



v(t) 




//// <& / 



X = 



\ \ \ 



+ C 



i(t) 



A. General Equations 



E a (t)i x 0<x<a 



E b (t)T x -b<x<0 



f E x d x =v(t) = E b (t)b + E a (t)a 



n . [j a - J b ] + ^ = => a a E a (t) - a b E b (t) + £ [s a E a (t) - s b E b (t)] = 



dt 



dt 1 



E b (t) = ^-E a (t)^ 



° a E a (t)-a b 



v(t) 



-E.(t)a 



dt 



s a E a (t)-s b ^-E a (t)a 



= 



' + ¥ldE f ^_ g bV(t) s b dv 
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B. Step Voltage: v (t) = v u (t) 





v(t) 


► 




t 



Then — = V 5(t) (an impulse) 
At t=0 

[ a b J dt b dt b K ' 
Integrate from t=0. to t=0+ 




E a (t) = 



a b V 



Ae 



-tA 



E a (t = 0) 



8 b V 



A = V 



e,b + e h a 



E.(t)= h PbV (1-e-^) 
a w a a b + a b a v ' 



e a b + e b a 



E„(t) = ^-E,(t)i 
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C. Sinusoidal Steady State: v (t) = Re [v e 
E a (t) = Re[E. e**] 
E b (t) = Re[i b e jfflt ] 



Conservation of Charge Interracial Boundary Condition 



a a E a (t) - a b E b (t) + — [s, E a (t) - s b E b (t)] = 



E a [cr a + jcoe a ]- E b [a b + JC0 8 b ]= 



E b b + E a a = V 



E h = 



V E a a 



b b b 



V E a a 



E a [a a + jroe a ]- t---£— [<? b + J»e b ]= 



, b b 

v J 



°a + j roS a + j^K + JC0 8 b ) 



V r -i 

= + jcoe b J = 



V 



j(OS b +a b jco8 a +a a [b(a a + j(D8 a ) + a(a b + jco8 b )] 



a su - 8 a E a - 8 b Eb 



(Sa^b "S b CT a )V 



[b(a a + j(0 8 a ) + a(a b + jra8 b )] 
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D. Equivalent Circuit (Electrode Area A) 



I = (a a + JC0 8 a )E a A = (o b + JC0 8 b )E b A 



V 



R a C a jco+l R b C b jco+l 



R. ^ — . R.. — 

3 a a A' b a b A 



a 



s h A 



C = " a ' ' r = — 
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